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ABSTRACT 

We consider the Mode Coupling Theory (MCT) of Glass transition for a Binary fluid. The 
Equations of Nonlinear Fluctuating Hydrodynamics are obtained with a proper choice of the 
slow variables corresponding to the conservation laws. The resulting model equations are 
solved in the long time limit to locate the dynamic transition. The transition point from 
our model is considerably higher than predicted in existing MCT models for binary systems. 
This is in agreement with what is seen in Computer Simulation of binary fluids. 



PACS number(s) : 64.70P, 05.60.+W, 64.60C, 47.35. +i,05.20.-y 



The phenomenon of glass transition has been studied widely in recent years using the 
mode coupling models for dense liquids. This involves a feed back to the transport properties 
of the liquid from the coupling of hydrodynamic modes in the liquid. The key aspect of the 
Mode Coupling Theory (MCT) of Glass transition takes into account the enhanced feed 
back effects to the viscosity due to the coupling of the slowly decaying density fluctuations. 
As the density increases beyond a critical value the system undergoes a dynamic transition 
to a non-ergodic state. The long time limit of the time correlation of density fluctuations 
is treated as an order parameter for structural relaxation. The transition to the nonergodic 
state is characterized by a nonzero value of the order parameter or the so-called Non Ergod- 
icity Parameter (NEP). For a binary mixture the self consistent MCT for Glass Transition 
and the dynamic instability resulting from the feedback mechanism has been studied by 
several authors 0-0 in the past. The model equations in these works predict the dynamic 
transition at a much lower density compared to the one where a structural arrest is seen in 
computer simulation of similar systems. This aspect of the Mode coupling model for the 
binary fluid was indicated in the computer simulations results reported in Ref. 0. In the 
present work we have constructed the MCT for the binary system from the equations of 
the Nonlinear Fluctuating Hydrodynamics(NFH) 0| with proper conservation laws. We 
consider the correction to the transport coefficient using the mode coupling approximation of 
dominant density fluctuations. This gives rise to the feed back mechanism on the transport 
properties. The possible dynamic transition allowed by the model equations is analyzed 
in terms of the solution of the resulting integral equations. Our main flndings are very 
different from the existing results 0. The theoretical prediction for the critical density of 
the dynamic transition now shifts to much higher density bringing it in better quantitative 
agreement with the computer simulation results. 

The slow variables for a binary fluid are partial densities Ps{x), s = 1,2 and total mo- 
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mentum density g{x). These can be defined microscopically as, 

Ps{x) = rrisY^ 6{x-R%(t)) 

2 Ns 

US) = EE^aW'^(^-^lW) (1) 



s=l a=l 



where R^i^) and Pio,{t) are respectively the position and i-th component of the momentum 
of the a-th particle in the s-th species. Ns is the number of particles in the s-th species in 
the mixture and the total number of particles is given by, N = J2s ^s- Generalized Langevin 
equation for the slow variables are 

^ + V.[%1+ 7..^ = ^s . = 1,2 (2) 

at p dps, 

-K7+^j + Ps^i^ + Lij— = fi (3) 

Ot p OPs OQj 

where the repeated indices are summed over - a convention that we follow in the rest of 
the paper, unless specifically mentioned. The '-fs/s correspond to the inter and Self-diffusion 
in the two species. We use the expression ^ 7^,3 = di'^^sisdj- On adding the equations (0) 
with s = 1, 2 for the two components, the continuity equation dp/dt + V.^ = for the total 
density p ( = pi + P2 ) is obtained. Thus the diffusion kernel in the first two equations of 
(0) for s = 1,2 should cancel. We will use the simple choice here that 721 = —711 = 7oV^ 
and 7i2 = —722 = 7oV^. Similarly the noise corresponding to density fluctuation satisfies 
61 + 62 = 0. Noise is assumed to be white Gaussian and related to the corresponding 
bare transport coefficients. The / and 6 are taken to be uncorrelated. The bare transport 
coefficient Lij can be expressed in terms of the shear and bulk viscosities rjo and Co in an 
isotropic fluid. The longitudinal viscosity is given by To = {4:rio/3 + Co)/ Po- 

The Free energy functional F that is used in constructing the equations of the Nonlinear 
Fluctuating Hydrodynamics plays an important role here. It has two parts, the kinetic 
and the "potential" respectively denoted by F^ and F^- In computing the kinetic energy 
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contribution, we follow the standard procedure due to Langer and Turski 0. 

Fk = lJ dx[g\x)/p{x)] (4) 

where the total density p = pi + P2 appears in the denominator of (^). For the potential 
part of the free energy there is an ideal gas contribution together with the interaction term. 
We use the standard form |^ in terms of the direct correlation functions 



1 f Ps{x\ 1 r 

Fu{p) = — / dxps{x) 1] - / dxdxfCss,{x - xf)Sps{x)Sps,{xf) (5) 

rris J Pos 2msmsi J 

From the dynamic equations (0) and we compute the time correlation function |Tl[ 
Cssi{x — xf,t — t') = < ips{x,t)4's/{xf,tf) > between the hydrodynamic variables ips and 
ipsh The normalized form of the density-density correlation functions, Rssi{q,z), is defined 
as Rssi{q,z) = Css/{q, z) / ^/XssXs/s/ where Xss/ is the equal time density correlation function 
between i-th and j-th species. The partial structure factor Sss' |]12| is defined from the 
Xss' through the relation Xss' = O'sO'sfSss/ where = msy/nl , nis and Us are respectively 
mass and particle density corresponding to the s-th species. We compute the mode coupling 
contribution to the transport coefficients by treating the non linearities in equations of motion 



with standard methods ||TT| . The renormalized form for the longitudinal viscosity at the one 
loop order is given by, 

1 r dk 

r(g,t) = — ^ / -—Vss,{q.k) Va,{q,h) QUk,t) CUki,t), (6) 

where ki = q — k. The vertex function Vssi{q, k) is given by, 

[msmsi\Vss,{q,k) = {q.k) Css,{k) + q.ki Css,{ki) (7) 

with Cssi{x — xi) denoting the equilibrium two particle direct correlation function |T^. The 
quantity Cssi{q) = {nsns,)^Cssr{q) is related to the structure factor Sssr{q) through the Ornstein 
Zernike relation [6is — Cis{q)]Sis,{q) = Sssr- In obtaining the above expression(^ we have 
considered only the coupling of the density fluctuations coming from the non linearity in 



the Pressure term that appears in the momentum conservation equation (^. This involve 
taking the density fluctuations as dominant and using the one loop correction or the so-called 
Kawasaki approximation to the four point functions to the transport coefficients. This is 
done in the same spirit of the self-consistent mode coupling approximation of taking density 
fluctuations to be dominant as in the case of one component fluid. We like to point out here 
that in the appropriate limit these results reduce to the one component fluid result that is 
in complete agreement with all other wave vector dependent models fl^ for one component 
systems. 

Next we consider the implications on the dynamic transition as a result of the new set of 
integral equations for the NEP. The ideal glass phase is characterized by the nonergodicity 
parameter, limj^oo Rss/il, t) = fss/il)- The correlation of density fluctuations is expressed in 
terms of the renormalized transport coefficients [0. In the asymptotic limit of long times we 
obtain the following set of self-consistent equations for the NEP's fssfi^), 

^-'^^^ - 1 + mm ■ 

Here T{q) is the long time limit of T{q, t) given by equation (|). VL{q) and 53(g) are given by, 

[Sss{q)Ss>si{q)\^ 

with Ass/ = CLsas,/po- Eq. (Q) constitute a set of coupled nonlinear integral equations for the 
non-ergodicity parameters fss/- The dynamic instability of the ideal glass transition in the 
binary system is then located from the self- consistent solution of eq. (|^) by iterative method 
in a similar manner as in the one component systems. The static structure factor goes as 
an input in the calculation. The interaction potential for the system enters the present 
theoretical description for the dynamics through this thermodynamic quantity. Due to the 
available simple formulation for the structure and in order to keep continuity with pervious 
works on binary fluids we will choose a mixture of hard spheres. This is described by 
three independent parameters, (a) The fractional concentration of particles Xs{= Ns/N), (b) 



The size ratio a{= cri/cr2), and (c) The total packing fraction rj = r]i + ri2 where rji and 1^2 
are the packing fractions of the individual species, ris = l^^sC^^, rig and as being respectively 
the number of particles per unit volume and diameter of the s-th species (s =1,2). We 
choose (72, diameter of bigger species, as the unit of length. The structure factors for the 
binary liquid required in computing the mode coupling vertex functions appearing in (H) are 
obtained using the standard results of the solution of Percus Yevick equations for the hard 



sphere ||15| , mixture. We solve for the nonlinear integral equations to search the nontrivial 
fixed points with non zero values of fssf{(l)- Depending on the thermodynamic parameters 
described above a nonergodic phase is seen beyond a critical density. The results remain 
unchanged if we replace x and a respectively by 1 — x and 1/a meaning an interchange 
of the label for the two species. In Fig. 1, we show the non- zero solutions for the NEPs 
fii and f22 respectively for the critical packing fraction 77^ = .59. The nature of the ideal 
glass instability in the binary mixture is dependent on the size ratio a as well as the relative 
abundance x. We list below a few comments. 

(i) The NEP f22 reach a small value at an intermediate wave number qm depending on x 
(fraction of the bigger particles 2). As the relative proportion x is increased this wave 
number become 0, beyond a value x > Xq. This behavior is primarily due to the 
nature of the structure function 53ss/(g) that appears in the mode coupling equation 
(|^) for the NEPs. In the inset of Fig. 1, we show how the position of q^ shifts with x 
at a = .7 and becomes zero at Xq = .977. 

(ii) The transition to the non-ergodic phase, indicating a structural arrest, shifts to higher 
densities as the size ratio (a) is decreased. For example, at x = .1 and a = .2 structural 
arrests are absent np to very high densities (77 < .64). This absence of transition here 
can be explained from the fact that for low values of x and size ratio a the smaller 
particles can diffuse easily through the cage of the bigger particles and hence the 
structural arrest is avoided. Indeed it is a competition between the two quantities. 



the size ratio a allowing easier movement of the particles and x giving the relative 
abundance of the voids formed between bigger particles. To stress the importance of 
the present work, it should be noted that, with the same values of the thermodynamic 
parameters the Ref. report a transition in the range (.51 < rj^ < .52). 

(iii) With the ratio of particle sizes, a smaller than .75, there is no structural arrest, over a 
range (.22 < x < .58). In the earlier mode coupling model equations this sensitivity 
to the sizes of the particles (a) and (x) was absent. The same qualitative dependence 
of a on the freezing in a binary fluid system has also been reported in studies related 
to the thermodynamics [|16| of the system. 

Fig. 2 shows the critical packing fraction {rjc) with x at constant a = 0.8. If the particles 
are similar to each other (a = 1) and x close to unity, the dynamic instability occurs at 
7] = 0.515 that is very close to the critical packing fraction ( with the use of Percus-Yevick 
Structure factor ) for a single component system as is expected to be the case. 

As a result of the mode coupling instability, the transverse current correlation function 
[§] also develops propagating shear modes [0. The expression for the speed of the shear 
waves (ct) in the binary fluid can be obtained in terms of nonergodicity parameters in the 
following form, 

= 7^\^ I dkk^^fu>ik)8fss>ik)Uik)Uk)JK^ (10) 
oOvr^^ y nin2 J 

where cfssf{q) denotes derivative of Css/{q) with respect to q. In Fig. 3, we show the behavior 
of Ct with (x) for two different size ratios(a). 

In the present work we have obtained, from the equations of Nonlinear fluctuating hy- 
drodynamics, the feedback mechanism of the self consistent MCT. Solutions of our model 
equations indicate that the dynamic instability shifts to much higher densities than what is 
predicted in the earlier works [^. This also agrees with the trend seen in computer simulation 
results^. The kinetic energy term, eq. (^, of the Free energy functional that is used here 
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is crucial for the construction of the model equations. This is also essential for obtaining 
in the momentum conservation equation the convective or the so-called Navier-Stokes term 
reflecting Galilean invariance. In considering the dynamic nonlinearities that gives rise to the 
mode coupling instability we have taken here a driving Free energy functional of Gaussian 
form. Thus origin of the nonlinearity is dynamic as is the case for the one component fluid. 
The higher order terms in the Free energy functional will give rise to nonlinearities in the 
density equations and this will require for the renormalization of the diffusion coefficients 
(7o). The lowest order vertex then is proportional to the 70 and will come from non-Gaussian 
terms in the free energy functional involving higher order direct correlation functions. We 
end the paper with few comments, 

(i) The earlier versions of the MCT ^ for binary systems can be reproduced from our 
work, using the same Gaussian Free Energy functional Fuas considered here, if one 
(a) replaces the kinetic energy term with i^fe = / dx[g\/2pi + g\l1p2\, §[ and (fz 
denoting the momentum densities of the individual components in the binary system, 
both treated as two separately conserved variables and (b) ignore the inter diffusion 7 
as well, (a) will also violate the Galilean invariance in the hydrodynamic equations. 

(ii) Here we have considered the dynamics of density correlation functions only in the 
asymptotic limit so as to study the implications on the dynamic transition. The utility 
of the present approach is demonstrated through the better result for the location of 
the dynamic transition. To study the dynamics over different time scales the time 
dependent mode coupling equations for the various density correlation functions has 
to be considered. These equations involve the bare transport coefficients 7jj and L and 



will be reported elsewhere [18 



(iii) Finally, the complete picture of the dynamics will involve computing the coupling to 
the current correlations that restores the ergodicity in the longest time scale |jl8 . 
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The present work treats the dynamics for the two component system with the correct 
set of slow variables and shows quantitative agreement of MCT results with computer 
simulation studies and remedies the theories presented in earlier works. 
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Figure Captions 

Fig. 1 : nonergodicity parameters /ii(g*)(sohd line) and /22(?*) (broken hne) at rjc — .59, 
x=.8 and a — .7. q*—qa2 is plotted along x axis. Inset shows ^'^(see text) Vs x plot for 
a = .8 and r] = .59 

Fig. 2 : critical value of rjc is plotted with x for a = 0.8. Regions "a" and "b" represent 
liquid and glass phases respectively. 

Fig. 3 : speed of shear mode ct is plotted as a function of x for two values of q;=0.85 
(solid) and 0.90 (broken) at 77 = .57. 
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